Thermal neutron flux produced by EAS at various altitudes 
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The results of Monte-Carlo simulations of Extensive Air Shower are presented to show the difference of hadronic 
component content at various altitudes with the aim to choose an optimal altitude for the PRISMA-like experi- 
ment. CORSIKA program for EAS simulations with QGSJET and GHEISHA models was used to calculate the 
number of hadrons reaching the observational level inside a ring of 50 m radius around the EAS axis. Then the 
number of neutrons produced by the hadronic component was calculated using an empirical relationship between 
the two components. We have tested the results with the ProtoPRISMA array at sea level, and recorded neutrons 
are close to the simulation results. 
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it more powerful and informative. 



1. Introduction 

A novel type of Extensive Air Shower (EAS) 
array (the PRISMA project) has been proposed 
some years ago il^ to study cosmic ray spec- 
trum and mass composition in the "knee" re- 
gion. Existing experimental data in the knee re- 
gion contradict each other and new approaches to 
the so-called "knee problem" are needed to solve 
this complicated and old problem. The PRISMA 
project is based on an idea that the main EAS 
component - hadrons has to be measured first of 
all. Special detectors (en-detectors) have been 
developed for this purpose. It was proposed to 
combine the central part of the PRISMA array 
with the LHAASO project [2^ by introducing the 
en-detectors in the center of LHAASO array. The 
en-detectors will make the LHAASO array sensi- 
tive to the hadronic EAS component thus making 
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2. Calculations 

We used the CORSIKA program [3] (ver. 
6.900) for EAS simulations with QGSJET and 
GHEISHA models. Calculations were performed 
for two primaries: proton and iron and for two 
altitudes: near sea level (170 m a.s.l.) and 
high mountain 4300 m a.s.l. As a first step we 
made calculations for fixed primary energies from 
lOTeV through lOPeV and zenith angles - 45°. 
We present here results for the number of hadrons 
and produced by them thermal neutrons inside 
a ring of 50 m radius around the EAS axis as 
a function of primary energy. The distributions 
over these numbers are also obtained. 

As any other Monte-Carlo program for EAS 
simulation, CORSIKA can not process particles 
with very low energies. We used following cuts for 
particles: 50MeV for hadrons, 0.5GeV for muons, 
60keV for electrons and gammas. The mean num- 
ber of evaporation neutrons < n > produced by 
hadrons in 3-m layer of surrounding soil and/or 
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construction materials has been calculated using 
an empirical relationship between them and the 
parent hadron energy (in GeV): 

<n> w36-£:^-5^ (1) 

This relationship originated from secondary par- 
ticles production in hadronic interactions and was 
obtained taking into account experimental data 
[1] and atomic mass A dependence of neutron pro- 
duction: < n > ^A°-^[5]. Therefore, the total 
number of produced secondary neutrons should 
be summarized over all hadrons: 

< ntot >= ^36-£;,°-^6 (2) 

It is easy to see that due to slow dependence on 
hadron energy and due to slow change of hadron 
mean energy with primary energy, eg. 2 can be 
simplified to: 

< ntot > »^36- < Eh >"-'^ < Eh >"-'^ (3) 

The latter means that the total number of evap- 
oration neutrons produced in EAS should be 
more or less proportional to the number of high- 
energy hadrons reaching the observation level. 
The great bulk of these neutrons are thermal- 
ized later. Thus recording thermal neutrons by 
detectors spread over big enough area one could 
recover the number of hadrons in EAS. This idea 
is a basis of the PRISMA project. 

3. Simulation Results 

Mean numbers of produced neutrons and 
hadrons inside a ring of 50 m as a function of 
primary energy are shown in Fig. 1 for primary 
proton and iron and for near sea level. As one can 
see, all dependencies can be fitted with power law 
functions. All indices for protons are close to 1: 
~1.14 at sea level and ~1.05 at 4300 m. For iron 
the indeces are little bit steeper at these energies. 
Small difference between primary proton and iron 
at the highest energy makes us sure that recovery 
of primary energy would be more adequate in the 
knee region. One can also see that the mean num- 
ber of neutrons at sea level and at 4300 m a.s.l. 
differ by a factor of ~ 10. The latter means that 
at Yangbajing level the array threshold energy 
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Figure 1. Mean numbers of hadrons and sec- 
ondary evaporation neutrons inside ring of 50 m 
radius as a function of primary energy at sea level 
for proton and iron primaries and at 4300 m for 
primary proton. 



can be lower by a factor of 10^/^'^ w 4. There 
are shown in Fig. 2 the distributions of different 
EAS components numbers inside 50 m around the 
axis for 2 primary energies: 30TeV and IPeV. An 
interesting issue here is rather high mean number 
of produced neutrons in comparison with hadrons 
and muons. The latter means that thermal neu- 
trons could be a better primary energy estimator 
than muons or high energy hadrons. Figure 2 also 
presents that at high altitude, within the primary 
energy range 30TeV-10PeV, the secondaries have 
higher quantities and then lead to higher energy 
resolution than ones at sea level. It is necessary to 
test these simulation results at different altitudes. 

4. Test at sea level 

4.1. The ProtoPRISMA array 

We have made several tests with the PRISMA 
prototypes in Moscow. Our test results are con- 
sistent with the simulations. A prototype of 
the PRISMA project array (ProtoPRISMA) [5] 
has been developed and started running on a 



3 




Figure 2. Different EAS components size distributions for primary protons and primary energies equal 
to 30 TeV (left panel) and 1 PeV (right panel) for sea level (upper panel) and for 4300 m a. s.l. (bottom 
panel) 
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base of the NEVOD-DECOR detector at National 
Research Nuclear University MEPhI in Moscow 
(170 m a.s.L). It consists now of 32 inorganic 
scintillator en — detectors situated inside the ex- 
perimental building at a level of 4th floor just 
around the NEVOD water pool. The neutron 
recording efficiency of the en-detector is equal 
to €d ~ 20%. The detectors have a cylindri- 
cal shape with the scintillator area equal to 0.36 
m? . It is now our standard en-detector made 
on a base of commercial polyethylene (PE) wa- 
ter tank of 200 liters volume. The scintillator 
thin sheets are situated on its bottom and are 
viewed by a single 6" PMT (FEU-200). The scin- 
tillator compound ZnS(Ag)-f LiF enriched with 
^Li up to 90% is very effective scintillator for 
heavy particle detection. It produces 160000 pho- 
tons per a neutron capture through the reaction 
^Li(n,a)t -f 4.8 MeV. It allows us to collect more 
than 50 photoelectrons from PMT photo-cathode 
per n-capture. Due to thin scintillator layer (30 
mg/cm?), a single relativistic charged particle 
produces very small signal. But, in a case of EAS 
passage, correlated signals from many particles 
are summarized and can be measured. There- 
fore, the same detectors are used to measure two 
EAS components: hadronic (neutrons) and elec- 
tromagnetic. That is why we call them as en- 
detectors (electron- neutron). The array consists 
now of 2 clusters of 16 detectors each. The clus- 
ters have their own triggering system and work 
independently. Coincidence of 2 or more hit de- 
tectors within time of 2 /is produces a simple trig- 
ger and on-line program preprocesses data, selects 
EAS events and marks them in accordance with 
different mathematical trigger conditions. All 
pulses are integrated with a time constant of 1 
/iS and are digitized using 10 bits flash ADC with 
a step of 1 /xs. In a case of powerful enough EAS, 
full pulse waveform of 20 ms duration is stored. 
For energy deposit measurement in a wide range 
we use 2 signals from each detectors: from the last 
12th dynode and from an intermediate 7th dyn- 
ode as well. Due to this we have the energy range 
more than 10^. Delayed pulses from thermal neu- 
trons captures are measured in each detector in 
the time gate 0.1 - 20 ms. 

In order to test our simulations at sea level. 



we have processed the experimental data accumu- 
lated for 2 winter months with the aim to check 
eq.(l) normalization factor. It is obvious that ab- 
solute value of neutron yield depends on the ex- 
perimental conditions. In our case the array is 
not situated in open air, but inside the build- 
ing. Therefore, the result can be regarded as 
an estimation only. Detailed neutron yield mea- 
surements will be performed later when the next 
PRISMA prototype will be running in Yangba- 
jing area in real conditions. 

4.2. Experimental results 

Mean number of recorded neutrons as a func- 
tion of EAS size (iVg) is shown in fig. 3 in 
comparison with full-scale Monte-Carlo simula- 
tions made using CORSIKA (with the parame- 
ters mentioned above) and the ProtoPRISMA ar- 
ray simulations. For neutron production we use a 
formula similar to eq.(l) but with a normalization 
factor taking the our neutron recording efficiency 
into account. The efficiency can be estimated as 
follows: e « Erf • e^m ■ Sarea, whcrc £em ~ 0.04 is 
a probability for neutrons produced in soil (con- 
crete) to escape to air [7 and earea=0.028 is a 
ratio of total detectors area (13.5 m^) to the ar- 
ray area (484.5 ) (coverage) . Taking these ef- 
ficiencies into account one can transform eq.(2) 
to: 

< ntot > « C • e-36-£:^-5'^w C-0.0080-£;^-^^ (4) 

This relationship was used in the simulations 
of the experiment. The normalization factor 
(C) was varied to adjust it to the experimen- 
tal points. Points in fig. 3 were obtained with 
C=0.625. Agreement is rather well (note that 
due to limited statistics for the highest EAS en- 
ergies corresponding calculated points at the right 
tail are little bit lower). As one can see, the 
mean number of recorded neutrons in the Pro- 
toPRISMA experiment is ~ 1.6 times lower than 
expected one. The reason could be found in the 
experimental details: detectors are placed inside 
building (not in open air), the structure of the ar- 
ray is asymmetrical and the detector coverage in 
2 clusters are dilferent [6] and finally, an existence 
of a large water pool inside the array affects un- 
doubtedly on the neutrons yield making it lower. 
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Figure 3. Mean number of recorded neutrons as a 
function of EAS size for Moscow site. Histogram 
- experiment, points - simulations. 



5. Summary 

In the frames of Russian-Chinese coUabora- 
tion Monte- Carlo EAS simulations have been per- 
formed and absolute normalization to the exper- 
iment has been done. We obtained rather low 
yield of neutrons at the ProtoPRISMA location in 
Moscow due to specifics of the experiment. Nev- 
ertheless, even at such yield and even at sea level 
it is possible to obtain some new results using the 
novel type of EAS array. We are looking forward 
to make such experiment at high altitude in Yang- 
bajing region. It would make us possible to check 
the method at high altitude, to make calibration 
using the ARGO-YBJ facilities and to measure 
real neutron yield in the future experiment site. 
The present calculations show that the number 
of produced neutrons at Yangbajing level is by a 
factor of 10 higher than that at sea level in EAS 
of equal energy. 
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